The flat X-ray spectrum of the LINER NGC 1052 
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ABSTRACT 

We report on ROSAT and ASCA observations of the LINER NGC 1052, which is 
the first one where broad optical lines in polarized light have been observed. The 2- 
10 keV spectrum is very flat, with a reserved — 0.1. A model where a nuclear source is 
- partly or totally - obscured by a screen of matter with column density ~ 10 23 cm" 2 
is the most convincing explanation for the observed flatness. This agrees with the 
hypothesis that the LINERs are a population of low-luminosity AGN, to which the 
Seyfert unification scenario applies. The intrinsic spectral index is still rather flat 
(rinstrinsic — 1.0-1.4), as observed in a few type-2 Seyferts so far or predicted if the 
accretion occurs in an advection-dominated flow. 
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1 INTRODUCTION 

It is still a matter of debate whether Low Ionization Nuclear 
Emission Line Regions Galaxies (LINERs) are powered by 
intense starbursts, or by an active galactic nucleus (AGN), 
or both. In some well-studied LINERs the non-stellar nature 
is almost certain (e.g.: M81, M87, M104). However, for the 
majority of objects the existence or relative amount of the 
AGN contribution remains unknown. Several facts point in 
favor of the idea that most LINERs are low- luminosity AGN: 
the similarity in the host galaxy properties and magnitude 
distribution (Ho et al. 1997), the discovery of supermassive 
black holes in the central regions (Ho 1998), and of broad 
wings of permitted optical lines (Ho 1999). 

Studies in the X-rays allow in principle to probe the 
innermost regions of the galaxies up to a few Schwarzschild 
radii around the putative central black hole and to inves- 
tigate the high-energy processes in the surrounding envi- 
ronment. Actually, soft compact X-ray emission has been 
detected in a handful of LINERs by the ROSAT HRI (Wor- 
ral & Birinkshaw 1994; Koratkar et al. 1995; Fabbiano & 
Juda 1997), but the instrumental spatial resolution allowed 
to limit the extension of the emitting region no better than 
about a few hundred pc. The ASCA observation of LINERs 
revealed that a "canonical" simple model is a general good 
description of the observed spectra, which is composed by a 
hard power-law with photon index V ~ 1.7-1.8 and a soft 
excess of thermal origin, with typical temperatures ~1 keV 
(see Ptak et al. 1998 and references therein). Although the 
latter component closely resembles the non-thermal power- 
law commonly observed in Seyfert Is, the lack of iron K Q 



line emission (Terashima et al. 1998) argues against this 
identification. 

If the LINERs are indeed a low-luminosity and nearby 
population of Seyfert-like objects, one might expect that 
the Seyfert unification scenario (see Antonucci 1993 for a 
review) applies to them as well and that a population of 
LINERs exists, where the Broad Line Region is hidden by 
optically thick matter. Actually, only a fraction of the ob- 
served LINERs shows broad optical lines (Ho 1999) . A sub- 
stantial support to this hypothesis has come from the recent 
discovery of a broad H a line in the polarized light spectrum 
of the prototypical LINER NGC 1052 (Barth 1998), con- 
firming previous weaker claims (Ho et al. 1997). 

If the above interpretation is correct, in some LIN- 
ERs one might observe X-ray spectra that are analogous 
to the ones typically observed in Seyfert 2 galaxies: either 
a power-law spectrum is seen in transmission through a 
Nh ~ 10 21 " 24 cm" 2 ; or, if the absorbing matter is optically 
thick to Compton scattering (Nh ~ 10 24 cm" 2 ), only the 
reflection of the nuclear radiation is detected. The last case 
has been observed in a handful of so-called "Compton-thick" 
sources (Iwasawa et al. 1997; Matt et al. 1996a, 1996b, 1997; 
Maiolino et al. 1998), including the LINER NGC 6240 (Iwa- 
sawa & Comastri 1998). 

NGC 1052 is an elliptical galaxy (3'x2'.l) at a distance 
of about 28 Mpc (or z = 0.0049, De Vaucouleurs 1991), with 
m B = 11.8. It is a LINER (Heckman 1980; Sadler 1987), lo- 
cated in the group HG 44 (Huchra & Geller 1982; Garcia 
1993), whose barycenter lies about 8' North from NGC 1052. 
Optical and UV observations show an extended (20") emis- 
sion line region with shock heating emission spectrum (Fos- 
bury et al. 1978; 1981) and an IR excess, which is confined 
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lable 1. Observation log 



Instrument Date Exposure time (ks) 

HRI 04/08/94 ~3.4 

HRI 02/10/95 ~22 

PSPC 23/07/1993 ~14 

ASCA 12/08/1996 ~36 (SIS) 

~38 (GIS) 



to a central region smaller than 2" (~ 300 pc) and coinci- 
dent with the optical nucleus. This suggests that IR emission 
could arise from dust heated by the non-thermal source ob- 
served in optical and UV (Becklin et al, 1982). NGC 1052 
is a strong radio source and its radio spectrum peaks at the 
GHz frequencies. The spectrum is flatter than a— -0.9 above 
1 GHz and is variable on a timescales of months (Disney & 
Wall 1977). VLBI observations showed a compact radio core 
with a diameter of 0.001 arcsec (~ 0.14 pc) and a radio halo 
of about 21 arcsec (Fosbury et al. 1978). Recently very lu- 
minous water masers have been detected toward this source 
(Claussen et al. 1998), which are not aligned perpendicular 
to the radio jet, like in NGC 4258 and NGC 1068, but rather 
lie along the jet at a distance of 0.07 pc from the central en- 
gine. This implies that masers emission in NGC 1052 should 
be associated with radio jet or radio continuum rather than 
with a molecular region orbiting around the central engine 
(like in NGC 4258). 

In this Letter we report still unpublished results on ob- 
servations of NGC 1052 taken by the ROSAT and ASCA 
X-ray observatories. Statistical uncertainties are quoted 
at 90% level of confidence for one interesting parameter 
(A^ 2 = 2.71); energies are quoted in the source rest frame; 
Ho = 50 km s _1 Mpc -1 and qo = 0.5 are assumed, unless 
otherwise specified. 



2 DATA PREPARATION 

Data of NGC 1052 have been retrieved from the ASCA and 
ROSAT public archives as screened event file lists, which 
have been pre-processed and reduced according to standard 
criteria. The log of the observations presented in this Letter 
is reported in Table [l] 

The ROSAT/PSPC is a Position Sensitive Proportional 
Counter sensitive in the 0.1-2.0 keV band with a spatial res- 
olution of about 20" for on-axis sources. Scientific products 
have been extracted from a circular area of 1', which includes 
99.5% of the 0.2 keV source photons (and an even higher 
fraction at higher energies). The background spectrum has 
been extracted from an annulus with internal radius of 3' 
and external radius of 5' draw around the source. In the anal- 
ysis of PSPC spectra we eliminated PHA channels below 9 
and above 200 due to calibration uncertainties. ROSAT/HRI 
is a High Resolution Imager with a spatial resolution of 
about 10" and a limited (2-band) spectral response. ASCA 
payload (Tanaka et al. 1994; Makishima et al. 1996) is con- 
stituted by a pair of Charge Couple Devices (Solid Imaging 
Spectrometers, SIS0 and SIS1), sensitive in the 0.6-9 keV 
band, and a pair of gas scintillating proportional counters 
(Gas Imaging Spectrometers, GIS2 and GIS3), sensitive in 
the 0.7-10 keV band. Only grade 0,1,2 and 4 of the SIS data 
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Figure 1. Radial profile of the combined NGC 1052 HRI ob- 
servations {circles), superimposed to the instrumental PSF, as 
determined from the white dwarf HZ 43 profile (solid line). 
The instrumental background (field of view average count rate 
5.4 X 10 -3 arcsec -2 ) is accounted by a constant plateau in the 
PSF. Each data point has a signal-to-noise ratio 5 

have been used. Scientific products have been extracted from 
circular areas of 4', 3.25' and 6' around the apparent cen- 
troid of the source, for the SIS0, SIS1 and GIS, respectively, 
including more than 99% of the source photons. To remove 
any possible contamination from HG 44, background spec- 
tra have been extracted using the complement area of the 
same chip where the source lies and from an annulus sur- 
rounding the source extraction region and having its same 
area for the SIS and GIS, respectively. The following results 
are not substantially affected if different regions of the detec- 
tors' field of views or spectra extracted from blank sky fields 
are employed. NGC 1052 is well above the background up 
to 8 and 10 keV in the SIS and the GIS, respectively. Spec- 
tra have been rebinned in order to have at least 20 counts 
per channel, in order to ensure the applicability of the \ 2 
test. Total count rates are: 0.0519 ± 0.0014, 0.0376 ± 0.0012, 
0.0403±0.0013, 0.0436 ±0.0015, and 0.0311 ±0.0017 s" 1 for 
the SIS0, SIS1, GIS2, GIS3 and PSPC, respectively. 



3 SPATIAL AND TIMING ANALYSIS 

In Fig. [I] the NGC 1052 radial profile extracted from 
the combined image of the two HRI observations is com- 
pared with the instrumental Point Spread Function (PSF). 
NGC 1052 shows extension up to 100", which is likely 
to be associated with the diffuse emission of the galaxy. 
The count rate associated with the extended component is 
0.168 s~ . It correspond to a 0.1-2 keV The luminosity of 
~ 3x 10 40 erg s~\ if one assumes a thermal bremsstrahlung 
with temperature 0.5 keV and seen through an absorbing 
column density of 3 x 10 20 cm -2 . No appreciable variability 
is detected, either in the ASCA 0.5-1.5 and 1.5-10 keV or in 
the PSPC 0.5-2.4 keV light curves. We will focus thereafter 
on the average spectral properties only. 
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1 2 
Energy (keV) 

Figure 2. Spectra (upper panel) and residuals in units of stan- 
dard deviations (lower panel) when a simple power-law model 
with photoelectric absorption is applied. Each data point has a 
signal-to-noise ratio > 5. Energies in abscissa are in the observer's 
frame 



4 SPECTRAL ANALYSIS 

Fitting separately the same model on the spectra of the 
ASCA detectors yields consistent results within the sta- 
tistical uncertainties. Particular care has been adopted in 
checking the consistency between the PSPC and the ASCA 
spectra. Several authors have observed a systematic differ- 
ence in the slopes measured by ROSAT and ASCA in the 
overlapping energy band, which ranges between 0.2 and 0.4 
(cf. Fiore et al. 1994; Iwasawa et al. 1998). Moreover, the 
source could have varied between the two observations, both 
in flux and in spectral shape, and the contribution of any 
diffuse emission from the underlying galaxy could be differ- 
ently spread underneath the source spot by the instrumental 
PSF. If we perform a simple power-law fit with photoelec- 
tric absorption in the 0.8-2 keV band, the spectral indices 
turn out to be consistent within the statistical uncertain- 
ties (Frosat = 1.5±q;5, Tasca — 1.1 ±0.3). Consistent re- 
sults with those presented in this Letter are obtained above 
0.5 keV if the ASCA data alone are fitted. At the light of the 
above considerations, we have fitted the spectra of all detec- 
tors simultaneously, allowing relative normalization factors 
between different detector types as free parameters in all the 
following fits. Their typical values are always very close to 
1 (N G is/Nsis = 1.02-1.07, Npspc/Nsis = 0.86-0.92). 

In Figure ^ the spectra and residuals are shown when 
a simple power-law model with photoelectric absorption by 
cold matter is applied. The fit is rather bad (x? = 1.59), 
mainly due to a prominent soft excess below ~1.5 keV and 
an apparent steepening of the spectrum above ~7 keV. A 
much better x 2 is obtained if a two-component model is 
adopted, where the soft excess is parameterized either via a 
thermal optically thick (blackbody) or thin (mekal) plasma 
or through a broken power-law. In all cases, an emission 
line is significantly required by the data (Ax 2 = 37), with 
centroid energy E c ~ 6.35 keV, consistent with K Q fluores- 
cence from neutral or mildly ionized iron. The 90% upper 
limit on the Equivalent Width (EW) of a fluorescent He- 
like (H-like) iron line is 170 eV (140 eV). Table § reports 



scribed in this section, which yield an acceptable x ■ All the 
phenomenological descriptions are characterized by a very 
flat hard spectral index (F ~ 0.1), which is totally incon- 
sistent with that observed in "canonical" spectra of LIN- 
ERs (Ptak et al. 1998). The observed fluxes in the 0.5-2, 
0.5-4.5 and 2-10 keV bands are ~ 3.2 x 10~ 13 , 1.3 and 
3.5 x 10~ 12 erg cm -2 s _1 , respectively, corresponding to un- 
absorbed luminosities of 3.6 x 10 40 , 1.3 and 4.7xl0 41 ergs -1 , 
respectively. It is interesting to note that the spectrum can- 
not be modeled with the composition of thermal components 
alone, regardless of their nature and/or different types of 
absorption. We have tried three different possible scenarios 
(albeit not mutually exclusive on physical grounds) to ac- 
count for the observed flatness: (a) a "partial covering" of a 
nuclear source; (b) a "pure Compton-reflection" model; (c) 
a "transmission" model, where the nuclear radiation is seen 
through a thick screen of absorbing matter. In cases (b) and 
(c) a soft excess is still required and has been parameter- 
ized with either a blackbody, an optically thin plasma with 
free abundance or a scattering component, i.e. a locally un- 
absorbed power-law with spectral index constrained to be 
equal to that of the primary power-law. A K a fluorescent 
iron emission line is required by all models as well. 

The "partial covering" yields a very successful fit 
(Xr = 0.95). A screen with N H ~ 10 23 cm" 2 , covering ~ 80% 
of the primary source is required. The photon index of the 
primary continuum is still very flat (~ 1.25). 

The "pure reflection" scenario is also a very good pa- 
rameterization of the spectrum (x? =0.90-0.93). The intrin- 
sic spectral index turns out to be even steeper than typically 
observed in Seyferts if the soft excess is modeled through a 
thermal component (r ~ 2.3), but perfectly consistent with 
it if the soft excess is due to scattering of the primary radi- 
ation. However, the iron line is rather tiny in this scenario 
(EW ~ 80-90 eV). 

The "transmission" scenario is again a good model for 
the data (x? = 0.95-0.97). If the soft excess is modeled 
through scattering or an optically thin plasma, the intrin- 
sic spectral index is again rather flat (F ~1.0-1.4) and seen 
through a Nh ~ 10 23 cm -2 screen. The unabsorbed lumi- 
nosity of the nuclear component is ~ 5.6 x 10 41 erg s _1 . The 
iron line EW ranges between 100 and 400 eV. 

We have also tried a "mixed transmission+reflection" 
scenario, assuming that the nuclear primary continuum 
is not a simple power-law, but is modified by Compton- 
reflection, probably originating in the putative relativistic 
Keplerian accretion disk spiraling around the supermassive 
black hole. In this hypothesis the iron line is expected to 
be significantly broadened (Matt et al. 1992; Tanaka et al. 
1995; Nandra et al. 1997), but the available statistics does 
not allow us to set any constraint on the line profile. An 
intrinsic Seyfert-like spectral index is obtained (F ~ 2), ab- 
sorbed through a screen of several 10 22 cm" 2 . However, 
this is got at the cost of an unplausible high relative nor- 
malization between the reflected and primary components 
(R > 20). The iron line is again rather weak (EW ~ 100 eV). 
It is very difficult to imagine a geometry such as to pro- 
duce this huge amount of reflection when the primary radi- 
ation is not completely suppressed in the ASCA bandpass 
by a Compton-thick absorber, unless the nuclear emission is 
strongly beamed towards the reflecting matter. 
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Table 2. Best fit parameters and x f° r the combined PSPC/ASCA spectral fits. Explanations of the model abbreviations: wa = pho- 
toelectric absorption; pc = partial covering; po = power-law; bp = broken power-law bb = blackbody; mk = MEKAL code (optically 
thin plasma); ga = Gaussian emission line; px = Compton-reflection 



Power-law + soft excess 



Model 

wa* (po+bb+ga) 
wa*(po+mk+ga) 
wa* (bp+ga) 



N H 
(10 20 cm- 2 ) 
o Q-l-l-6 
z - ci=t 0.8 

3±2 
4 1 + 1 - 6 



hard 



0.11± 



0.11± 



1.0. 08 
c 0.07 

0.15 ±0.07 



kT or E break 
(eV) 
165±J| 
490 ± 150 

1 40J_0.16 



r soft or z 



o.o7±;; 

1.9±8 



E c 
(keV) 

fi or i 0.07 

a « I 0.07 
D.OOICq og 



EW 

(eV) 
220±1J° 
240 ± 80 
240 ± 80 



X 2 / dof 

313.3/310 
310.1/309 
313.4/310 



Model 

wa*pc(po+ga) 



(10 



N 1 
20 c H m - 2) 



2 6+ ' 9 



1.21 ±0.16 



Partial covering 
N| C f E c 

(10 22 cm" 2 ) (keV) 

10.3±};| 0.79 ±0.05 6.36 ± 0.06 



EW 

(eV) 
230±« 



X 2 / dof 
296.9/310 



Model 

wa*(px+po+ga) 
wa* (px+bb+ga) 
wa* (px+mk+ga) 



(10 2 



5.0±° 
3.5±i 
4.3±J 



o 1R-I-0.16 
Z - 1D=C 0.14 
2.43 ± 0.09 
2 42+0.10 



Pure rcflcction+soft excess 



kT 

(eV) 

190 ± 30 
620±}™ 



Z or f s X (2w/n) 



n 0033+ 0011 
u.uu^di,, 0016 



0.05±g;g2 



E c 
(keV) 
6.4+ 
6.4+ 
6.4t 



EW 
(eV) 
80± 3 ™ 
90 ±70 
90±"!° 



X 2 / dof 

290.8/312 
283.0/311 
278.0/309 



Model 

wa* [wa* (po+ga) +po] 
wa* [wa* (po+ga) +bb] 
wa* [wa* (po+ga) +mk] 



(10 20 



N 1 



2 6± ' 9 
< 1.8 

2 X.+0.8 



1.22 ±0.16 
0.49±°; 2 ° 

i.o±8-5 



Transmission + soft excess 



N 2 

(10 22 cm- 2 ) 
12 ±2 

9±« 



kT 



260 ± 50 

> 9 



f s or Z 
(keV) 

U - ZD=C 0.09 



6.36 ±0.06 
6.36±°;jg 

fi i 0.07 



EW 
(keV) 
300 ± 110 
260 ± 90 
220 ± 110 



X 2 / dof 
(eV) 
297.0/310 
299.3/309 
296.0/309 



Model 

wa* [wa* (po+ga+px) +po] 



N 1 

(10 2() eV 2 ) 



5-5±l; 6 7 



2 1± ' 2 



Mixed transmission + reflection 
fs 



N 2 

(10 22 cV 2 ) 



3.0±i;; 



0.5±° 



R 

80±? 



E c 
(keV) 
6.38±H1 



EW 

(eV) 

ioo±l° 



X 2 / dof 
280.6/309 



t fixed 



5 DISCUSSION AND CONCLUSIONS 

The X-ray spectrum of the LINER NGC 1052 is ex- 
ceptionally flat. The observed 2-10 keV photon index is 
r ~ 0.1. The 2-10 keV luminosity in the same energy range 
(~ 5 x 10 41 erg s _1 ) is more than one order of magnitude 
higher than the expected contribution from unresolved Low 
Mass X-Ray Binaries within the galaxy (Canizares et al 
1987; Matsushita et al. 1994). Models to explain such a 
spectrum fall in two main classes: those which require a 
highly absorbed (Nh — 10 23 cm~ 2 ) transmitted power-law; 
( "partial covering" and "transmission" scenarios) ; and those 
where the 2-10 keV flux is dominated by Compton-reflection 
( "pure reflection" scenario). 

In the former cases, the spectral index is still rather 
flat (T ~ 1.0-1.4). Such values are distant from the average 
observed in Seyferts (Nandra et al. 1997) , although not un- 
common among type-2 objects (Smith & Done 1996; Turner 
et al. 1997). Alternatively, such a flat index emerges nat- 
urally if accretion occurs via an advection-dominated flow 
(ADAF, Narayan & Yi 1995). ADAF has been recently in- 
voked to explain the emission from low- luminosity AGN, 
thanks to the very low radiative conversion efficiency that 
characterizes it. It has been rather successful in modeling 
the Spectral Energy Distributions of some low-luminosity 
AGN (see e.g.: Narayan et al. 1998). It has been recently 
pointed out that the canonical model significantly overesti- 
mates the radio and sub-millimeter flux in a number of ellip- 



tical galaxies, hence suggesting the existence of mechanisms 
to suppress the emission from the central region at these 
wavelengths (Di Matteo et al. 1998). In the basic ADAF 
scenario the X/7-ray emission should be dominated by ther- 
mal bremsstrahlung due to the very hot electron distribution 
in the disk, with typical kT ~ 100 keV. In low-resolution, 
band-limited detectors such a spectrum might be well ap- 
proximated by a simple power-law with F ~ 1.4. Further 
observations covering the energy band up to a few hundreds 
keV might be very valuable is assessing the viability of such 
a scenario, and are nowadays possible thanks to the un- 
precedented sensitivity of the Phoswitch Detector System 
on board BeppoSAX. 

If the transmission scenario is true in NGC 1052 (either 
in the Sakura-Sunayev or in the ADAF flavor), two more in- 
gredients are needed: i) a screen of cold absorbing matter of 
Nh — 10 23 cm- 2 , which could be provided by the "torus" en- 
visaged by the unification scenario or (in the ADAF case) by 
the cold outermost regions of a geometrically bloated disk. 
The observed iron line is consistent with being produced by 
the same absorbing screen, even if a contribution to the iron 
line comes from a relativistic, X-ray illuminated disk close to 
the nuclear black hole (note that an ADAF should not intrin- 
sically produce any fluorescent line); ii) a soft component, 
either as partial covering, or as a thermal emission peaking 
at several keV or as an energy-independent ~ 25% scatter- 
ing of the primary radiation. The discovery of broad optical 
polarized lines advocates for scattering indeed to occur in 
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X-ray component from the HRI data, which is likely to be as- 
sociated with the diffuse emission of the galaxy. The limited 
spectral resolution of the HRI does not allow a self-consistent 
determination of the luminosity of such a component. If one 
assumes that it is thermal in origin with a temperature of 
0.5 keV (comparable to those typically observed in elliptical 
galaxies; cf. Fabbiano 1989), the implied 0.1-2 keV lumi- 
nosity is ~ 3 x 10 40 erg s _1 . From the observed correlation 
between the blue and the X-ray luminosities in ellipticals 
(David et al. 1992), one may expect a contribution to the 
0.5-4.5 keV luminosity from the diffuse galaxy emission of 
the same amount (David et al. 1992). Of course different 
components could be simultaneously present in the soft X- 
ray band, thus confusing the interpretation of the data. A 
detailed modeling of the soft excess in NGC 1052 requires 
experiments with much higher spatial resolution and/or ef- 
fective area. 

Alternatively, a Compton reflection-dominated spec- 
trum from a "canonical" Seyfert-like primary spectrum pro- 
vides an equally good parameterization of the data. How- 
ever, in this case the iron line is very faint, the EW being one 
order of magnitude lower than observed in Compton-thick 
Seyfert 2s so far (Iwasawa et al. 1997; Matt et al. 1996b; 
Maiolino et al. 1998) and expected on theoretical grounds 
(Ghisellini et al. 1994; Matt et al. 1996a). In principle, a 
strong iron under abundance could resolve the discrepancy, 
although it is not clear why this should be the case in the 
nucleus of an elliptical galaxy. Moreover, there is currently 
no observational evidence for iron under-abundance in AGN 
central regions. A more appealing solution may be that the 
surface of the inner disk is ionized of such an amount that 
the bulk of the line flux is suppressed via resonant trapping 
(Matt et al. 1993). We would therefore be observing only a 
minority of the line flux, emerging from the outermost and 
cold parts of the disk. However, it is difficult to imagine a ge- 
ometry where the inner regions of the disk are visible, while 
the primary nuclear emission is totally suppressed. 

While a wide literature exists on the optical spectra 
of LINERs, the X-ray spectrum of only a few objects has 
been studied in detail. In this Letter we have presented evi- 
dence that an obscured (Nh ~ 10 23 cm~ 2 ) AGN, possibly 
accreting via an ADAF, is the ultimate source powering 
NGC 1052. Undoubtedly, future X-ray observations will con- 
tribute to further enlighten the connection between LINER 
phenomenon and nuclear activity, and the nature of the ac- 
cretion processes occurring in LINERs. 
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